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ABSTRACT 


A review of systematic perturbation procedures for the analysis of non- 
linear problems is presented. The cases when the multiplicity of an eigen- 
value is finite or infinite are treated for self-sustained and for forced oscil- 
lations. The possibility of the formation of shock waves is discussed. Appli- 
cations to acoustic problems are presented. 


INTRODUCTION 


Most of the problems in acoustics can be treated successfully by the lin- 
earized theory since the nonlinear terms in the governing equations are in gen- 
eral of higher order. The linearized theory predicts the natural frequencies. 
When there is only one eigenfunction associated with a natural frequency, we 
say the eigenvalue is simple and the mode of the free oscillation is given by 
the eigenfunction times a constant, the amplitude a, which is undefined. For 
forced oscillations the amplitude becomes infinite when the forcing term is in 
resonance. 

When the small nonlinear terms are included, the periodic solutions can 
be constructed by the perturbation theory. The linear problem will yield the 
first term in the perturbation expansion of the solution, and further terms will 
also be determined by linear problems. This theory is based upon the original 
discovery by Lindstedt and Poincare that, to avoid the occurrence of secular terms 
in applying perturbation theory to periodic motions in celestial mechanics, it 
is necessary to make a perturbation expansion of the period of motion. This 
theory is frequently used for problems involving nonlinear ordinary differential 
equations. It has been applied systematically to periodic free vibration prob- 
lems involving nonlinear partial differential equations by Kdller and Ting 
(ref. 1). The solution and the frequency are assumed to be regular functions 
of a small parameter e, the order of magnitude of the nonlinear terms. 

When the eigenvalue is simple in the linear problem, the orthogdnality of 
the eigenfunction with the inhomogeneous terms in the governing equation for 
the next order solution, and those for the higher ones, removes the secular terms 
and also defines the amplitude of the oscillation as a function of the fre- 
quency. 
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When the multiplicity of an eigenvalue in the linear problem is equal to 
k, there are k linearly independent eigenfunctions. The linear solution is a 
linear combination of those k eigenfunctions. The orthogonality of the inhomo- 
geneous terms of the governing equations for the higher solutions with k eigen- 
functions yields k equations to relate the k coefficients in the linear combina- 
tions to the frequency. The nonlinear analyses define not only the amplitude 
but also specify the mode of the linear oscillation to one special linear com- 
bination of the k eigenfunctions. 

When the multiplicity of an eigenvalue is infinite, the linear solution 
can be represented as an infinite series of the eigenfunctions. The orthogona- 
lity conditions with the eigenfunctions yield infinite numbers of equations for 
the coefficients in the series. The solution of these equations becomes very 
difficult, in general. Of course, we can construct an approximate solution 
by truncating the infinite series to a finite number of terms and imposing only a 
finite number of orthogonality conditions. This approximate solution will be 
useful only when the infinite series happens to converge very fast. 

Fora slightly nonlinear wave equation without dispersion, the multiplicity 
of each eigenvalue is infinite. The general solution of the linear wave equa- 
tion can be represented by forward and backward waves with unknown wave forms. 

The orthogonality conditions were shown to be equivalent to an integral con- 
dition by Keller and Ting (ref. 1) and by Hale (ref. 2). Using the integral 
condition Fink, Hall, and Khalili (ref. 3) showed that it leads to a functional 
differential equation for the wave form and obtained explicit solutions in 
terms of elliptic functions for three types of nonlinear forcing terms. 
Generalization of the analysis and the establishment of an integral solvability 
condition for n-dimensional space have been made by Ting (ref. 4). 

The systematic perturbation theory was applied to several interesting non-, 
linear boundary value problems by Millman and Keller (ref. 5). They also pre- 
sented a systematic procedure for the construction of solutions of forced oscil- 
lations. The amplitude of the forcing term and the energy of the system are 
assigned. The solutions and the frequency are again regular functions of the 
small parameter e . The energy equation guarantees the boundedness of the solu- 
tion. With a simple eigenvalue, we obtain the finite amplitude solution of free 
oscillation when the forcing function is in resonance while its amplitude van- 
ishes. When the eigenvalue is not simple, the force function can be any linear 
combination of the eigenfunctions with assigned amplitudes, and the solution of 
the forced oscillation will remain finite due to the energy condition. However, 
in order to recover the finite amplitude solution of free oscillation, the forc- 
ing function has to be a linear combination of all the eigenfunctions with all 
the coefficients approaching zero simultaneously while their ratios, which 
remain constant, are specified by the energy equation so that the wave form of 
the forcing function will be related to that of the free oscillation. 

The procedures for the construction of the perturbation solutions and the 
statements regarding the solutions will be demonstrated in the following two 
sections for slightly nonlinear one dimensional wave equations with or without 
a finite dispersion term. When there is a finite dispersion term, the eigen- 
value of the linear equation is simple. Without the dispersion term, the 
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multiplicity of the eigenvalue is infinite. For each case, we will treat the 
problem of the free oscillations first and then that of forced oscillations. 

The author wishes to acknowledge Prof. Joseph B. Keller for valuable dis- 
cussions regarding this paper. 

SYMBOLS 

A amplitude of the forcing function 

a amplitude of the oscillation 

b. coefficients of linear combinations in eq. (31) 

vi 

coefficients of in eq. (15) 

E energy of the system 

f nonlinear term 

g, 6 the wave form of the solution and the forcing term respectively 

j, k, ni, n positive integers 
t time 

u the solution of the nonlinear equation 

X space variable 

a dispersion coefficient 

3 defined in eq. (28) 

7 coefficient in the nonlinear term when f(u) = yu^ 

e small dimensionless parameter denoting the order of magnitude 

of the nonlinear terms 

X eigenvalue or co^ 

5 the phase variable 

T integration variable 

u frequency 

(•) differentiation with respect to e at e = 0 

Subscripts: 

0 the leading term, i.e. , when s = 0 

x,t partial derivatives with respect to x or t 
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NONLINEAR WAVE EQUATION FOR A DISPERSIVE SYSTEM 
Analysis of Free Oscillations 


Let us construct the periodic solution of the nonlinear equation with a 
finite dispersion coefficient, 

Utt ■ ^xx ^ " ef(u) 0<X<1T (1) 

subjected to the periodicity and boundary conditions, 

u(x, t + 2Tr/(o) = u(x,t) (2) 

u(0,t) = u(w,t) - 0 (3) 

e is a prescribed small parameter and e f(u) represents the nonlinear force. 

0 ) is an undetermined angular frequency. Let us introduce t'= ut and u"(x,t') 


= u(x,t) in eqs. (1X2K3) so that the period in the new time variable is 2it. We 
will drop the primes and eqs. {1)0, 0) become 

w^u^^-u^^+au = ef(u) (4) 

u(x,t + 2ir) = u(x,t) (5) 

u(0,t) = u(ir,t) = 0 (6) 

We shall seek a solution u(x,t,e) and a corresponsing angular frequency w (e) 
which are representable by finite Taylor series in e, i.e. 

u(x,t,e) = u^(x,t) + eu(x,t) + ^e^u +... (7) 

and x(e) =co^(e) X +.... (8) 


where Uq = u(x,t,o), to = o)(0) and (’) denotes differentiation with respect toe 
at e = 0. ° 


By setting £ = 0, 
solution L^, i.e.. 


to^^ U„, . 

O ott 


eq. (4) becomes the linear equation for the zero-order 


- u. 


oxx 


a u„ = 0 
0 


(9) 


and the same periodicity and boundary conditions eqs. (5), (6) hold for u . By 
choosing the origin of t appropriately we can require u. = 0 at t = 0 an8 write 
these solutions as ^ 


u^ = a sin nx cos t (10) 

Xo = u^o =a+n2,n =l, 2, (11) 


For each integer n, or the eigenvalue x^ is defined. When a is an irra- 
tional number we can snow that there is only%ne eigenfunction, namely sin nx 
cos t. The amplitude a is undetermined so far. 
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We differentiate eqs. (4), (S, (6) with respect to e and set e = 0 to obtain 
the equations for the next order solution u. They are: 

and eqs. (5) and (6) with u replaced by u. They will have a solution (ref. 1) 
if the inhomogeneous part is orthogonal to the eigenfunction, i.e. 

[-Iuq + fCUp)! sin nx cos t = 0. (13) 

tt 


We will give a physical meaning to this condition. The inhomogeneous term in the 
linear equation for u can be considered as the forcing term. The solution will 
be finite only when the forcing term is not in resonance with the normal mode of 
the homogeneous system. In other words, the coefficient of the Fourier compo- 
nent, sin nx cos t, of the forcing terra should vanish as expressed by condition 
(13). 



The solution is 

OO 00 

U = S Z' C.. sin j X cos kt/[j^ + a - 
k=0 j =1 

The prime over the summation sign means j ^ n when 
are defined by: 


00 CO 

f(Uj^) = 1:1: 
k=0 j=l 



cos kt sin Jx 


(14) 

k = 1. The coefficients 
0<X<TT (15) 


The denominator in eq. (14) can be rewritten as - n^k^ + a(l-k^) which will 
not vanish since a is an irrational number and j ^ n when k - 1. Eq. (13) in 
turn yields the first order frequency-amplitude relationship, 

( 1 )^- eX+ O(e^) = - Cp /a + OCe^) (16) 

For the special case of f(u)= - u^, we have 

= -(9/16)a3 

and 0 )^ =«+ n^ + e(9/16)a^ +0(e^) (17) 

Figure 1 shows the amplitude-frequency curves fora '/T, s = 0.1 in the neigh- 
borhood of the first two natural frequencies, i.e., D = 1,2, 
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Analysis of Forced Oscillations 


We will modify the preceding problem by adding the distributed forcing 
function A sin j x coswt where j is a positive integer and A is the amplitude 
of the forcing function. The differential equation, eq. (4), must be changed 
to: 

M^u - u +«u - ef(u) = A sin jx cos t (18) 

u Xa 

while eqs. (5) and (6) remain unchanged. 

We will now establish the energy equation. We multiply eq. (18) by u., 
integrate with respect to x from 0 to ir and carry out integration by parts ^ 
for the second term with the aid of eqs. (5) and (6). We obtain 



u| + u^ + u2a 
t 


u 

2e J f(u)du 
0 


- A u(k,T) sin jx cos t dt}dx = 0 
0 

Therefore the energy E is a constant and in particular at t = 0, we have 

h f ^ *^x ^ j dujdx = E 

0 0 't=o 


(19) 


We will prescribe the energy E and consider u and u to be functions of E, A 
and e . We will then represent u and to as finite Taylor series in e as in 
eqs. (7) and (8). Letting e = 0 in eqs. (18), (19), (5) and (6), we obtain 


w^^u„+. - u^„„ +au^ = A sin jx cos t 
0 ott oxx 0 


and 


H. 




/ 

0 


+ u^ + 
“t “ 


.u^)| 


dx = E 

t=0 


( 20 ) 

( 21 ) 


while eqs. (5) and (6) hold for u^. The solution of eqs. (20), (5) and (6) is 
u^ = a sin jx cos t (22) 

where a = A/[j2+ a -w 2 ] (23) 

w is related to the energy E by eq. (21), 

0 


E - ira^( j^ +®)/ 4 

= (ir/4)A^(j^ +“)/(j^ +tt-a)Q^)^ (24) 

Since is an irrational number, the forcing term is in resonance with the lin- 
ear system only at the jth natural frequency. With = j^ + “ , the ampli- 
tude a and the energy E in the linearized theory become infinite. For the non- 
linear theory, we prescribe the amplitude A and the energy E and define a)^A,E,e), 
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from eq. (24). We will then proceed to determine a>^( A, E,e) from the next order 
analysis. 


We differentiate eqs. (18) and (19) with respect to e and the set e = 0 to 

I 


{ “ox“x - { “ 


We make use of the fact that Uq+ = 0 at t = 0. Again eqs. (5) and (6) hold also 
for u. The solution of eqs. (Zo), (5), and (6) is 


+^il 

u = ^ - - sin j X cos t 


Sin k X cos mt 


m = 0 k - 1 k2 + <x 


The symbol (') means that k ^ j when m = 1. Here the coefficients are de- 
fined by eq. (15). Insertion of eq, (27) into eq. (26) yields 



where e = u^(x,0) = a sin jx. Use of eq. (28) In eq. (8) yields in terms of 
A and E to the first order in e. 

Fromeq. (24) it is clear that for a nontrivial solution, a / 0 and E is posi- 
tii/e. We can solve for to 

O 

0)2^ +<X ±[TTA2(j2 +«)/(4E)]^ (29) 

In the limit A-» 0, we obtain the results for free oscillations, 

j2 + 0( 
and 

X -»■ - C . i / a 

since f(u ) is finite. These results are in agreement with the results of 
eqs. (li),u4) and (16) obtained in the preceding subsection by a different proce- 
dure. 
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Since the amplitude a of the solution is equal to {4E/[iT(j2 +o()3}'^ we 
win plot a vs 0)2 for each constant forcing amplitude A instead of E vs 
Fig. 2 shows the curves for A = 0, 0.2, 1 with oc = e = 0.1, j = 1, f(u) 

= - u^ Of course, the curve for A = 0 agrees with the curve in Fig. 1 near 
the first natural frequency. 

We will now explain why these two apparently different methods for the con- 
struction of finite amplitude free oscillations are equivalent. 

In the first method, the underlying principle is that the solution u or 
Uqj u etc., have to be finite. The orthogonality condition eq. (13) guarantees 
that u is finite since the eigenvalue is simple and the solution it in eq.(14) 
is therefore not in resonance with the natural mode which is proportional to u . 
In the second method, we note that Uq is finite for finite energy E. Hence 
fCug) is finite, and the first order energy equation (26) says that the part of 
Ux which is orthogonal to Uq„ has to be finite. Uq is an eigen-solution. 
Since^the eigenvalue is simpTe, eq. (26) is sufficient to guarantee 
that ii is finite. When the eigenvalue is not simple, and the forcing function 
contains only one of the eigenfunetions, eq. (26) is not sufficient to i 
insure that u is finite although eq. (26) will produce a relationship between 
and E for A = 0. Of course, if we continue to the second order energy equa- 
tion, the appearance the term would then require u to be finite. The proce- 
dures for handling the problems with nonsimple eigenvalues are described in 
the next section. 


NONLINEAR WAVE EQUATION FOR A NONDISPERSIVE SYSTEM 


Analysis of Free Oscillation 


The governing equations (4), ($, (6) in the preceding section remain appli- 
cable when we set «= 0. The perturbation expansions eqs. (7) and (9) are the 
same. The governing equations for the zero order solution are 




0 ^ott 


- u 


oxx 


= 0 


(30) 


and eqs. (5) and (6) with u replaced by Uq. Again we choose the origin of t- 
such that u^ = 0 at t = 0, the solution for a given integer n can be written as 

00 

u = a (2 b. sin nj x cos njt] (31) 

° j=l ^ 

with Xq =w2q = n2 (32) 

Uq is expressed in terms of a linear combination of all the eigenfunctions 
associated withXg- Since we factor out the amplitude a, we can impose a con- 
dition on bj's, say 

2 b2. 4^ = 1 (33) 

1 J J 
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The governing equations for the next order solution U are eqs. (12), (5), and (6) 
withc< = 0. The solution u is bounded when the inhomogeneous part of eq. (12) 
is orthogonal to all the eigenfunctions sin njxcosnjt for all j. The results 
are 

P 2ir ir 

lab. + ^ / dt f dx f(u„) sin njx cos njt = 0 (34) 

J TT 0 t) ° 


for j = 1, 2. — . Eqs. (34) and (33) are the equations for all the b.‘s and a 
in terms of X. They are difficult to solve. ^ 

It was observed in ref. 4 that at least for the special case of f(u) 

= 7 u^ where 7 is a constant, the approximate solution which was obtained by 
keeping only two terms in eq. (31) and applying eq. (34) for j = 1 and 2 differs 
from the exact solution by less than 0.1%. The construction of the approxi- 
mate solution is the same as if the multiplicity of the eigenvalue were finite 
(say equal to 2). For a general nonlinear problem, it would be advisable to 
construct the exact solution by the following procedure. 


The general solution of eqs. (30), (5), and (6) with X - n^ 
written as ° 

u(x,t) = ^a [g(t + x) - g(t - x)] 

with 

g(5 + 2ir) = g(?) 


= 1, can be 

(35) 

(36) 


where g is the unknown wave form. We introduce an extra factor a so that a 
normalization condition on g can be introduced, say 

2ir 

i_ j [g"(5)l2dg = 1 (37) 

TT 0 

The series eq. (31) can be identified with g as 

00 

g(g) = 2 bj. Sin j?. (38) 

1 

The orthogonality condition of eq. (34) can now be replaced by an integral equation 
for the wave form g (ref. 1, 2) 

- g(t) +/ f[u (t-x, x)] dx = 0 (39) 

d 0 ° 

For the case of f =7u^ , eq. (39) is reduced to a functional differential 
equation for g and the solution is an elliptic function (ref. 3). Conse- 
quently all the coefficients b. in eqs. (38) or (31) are defined. Eq. (39) also 
provides a relationship between a and X, i.e. , the amplitude frequency rela- 
tionship. 


Analysis of Forced Oscillation 


If we introduce a distributed forcing term proportional to an eigenfunction 
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(40) 


say the first one sin x cos t, eq. (30) becomes 

u„^.. - = A sin X cos t 

0 ott oxx 

The solution of eqs. (40), (5), and (6) is 

u„ = a sin x cos t (41) 

0 

where a = A/[l -w (42) 

0 

The energy equation is eq. (21) with W = 0. Letting e = 0, we relate to E 

£ = Tra2/4 = (tt/ 4) A^/d- (43) 


0)2 ;= 1 ± [ttA2/ {4E)]^ 

When the amplitude of the forcing function A and the energy E are non- 
zero, the amplitude a and the frequency o>q are defined by eqs. (42) and (43). 

In particular Wq f 1; therefore, a' is finite. 

The next order differential equation and energy equation for u and X are 
given by eqs. (25) and (26) with oi- 0. Their solutions are given by eqs. (27) 

and (28) with « = o and j - 1. In particular for f(u) = - u^* we have 

X cos t- • y5rir9; ; j' '' 2T sin 3 x cos 3 t (44) 

^ 16(1-0) 2) ' sin 3 X cos t + le f fT a T^) ^ x cos 3 t 

and 

X = (9 a2/16) [1 + 3 o )^2 (1 _ (1 . (45) 

In the limit of A ^ 0 and o)o^l such that the, energy A and also the ampli- 
tude a remain constant, we obtain X" 9a^/16 while u becomes unbounded because 
of the last term in eq. (44) unless E = a = 0. 

In order to obtain nontrivial free oscillations as a limiting case of the 
forced oscillations, the forcing function should contain all the eigenfunctions 
in the form of an infinite series eq. (31) or be represented as general solu- 
tion eq. (35) for X =1. Eq. (40) should now be 

o)^Q % ■ %x ^ (A/2 )[G(t+ x) - G(t - x)] (46) 

^tt 

The solution u^ is 

Uq = (a/2 .)[g(t + x) - g(t - x)] (47) 

where a = A/(l-o)^2) (48) 
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g" ( 5 ) = - 6(1) 


(49) 


with g (5 + 2ir) = g (c) and g(0) = 0. Since we introduce an 

extra factor A we can normalize G or g by eq. (37). The energy is related to a 
or A by eq. (43). 

In the limiting case of -»■ 1 and A 0, while E and a remain finite, 
u in eq. (47) is of the same foPm as eq. (35) for the free oscillation. The 
governing equation for Ci is also the same as that for the free oscillation. 

On account of the next order energy equation, the part of u which is infinite 
has to be orthogonal to all the eigenfunctions. Therefore, u is finite and the 
wave form g(g) will be the same as that for the finite free oscillation in 
ref. 3. The wave form for forcing term is in turn defined by eq. (49) as A -> 0. 


CONCLUDING REMARKS 


The perturbation theories for nonlinear free and forced oscillations are 
reviewed. In the case that the eigenvalue for the linear problem is simple, the 
solution for the forced oscillation yields that of free oscillation when the 
forcing term is in resonance while its amplitude approaches zero accordingly. 
This statement remains true when the multiplicity of the eigenvalue is finite or 
infinite provided that the forcing function is a linear combination of all the 
eigenfunctions. Their amplitudes vanish in such a menner that the wave form of 
the forcing function is defined by that ft)r the free oscillation. 

From eq. (24) we see that for a special combination of E and A, k^+a- 
vanishes for a pair of (k,m) and the solution u given by eq. (27) becomes in-° 
finite. We should then add a free oscillation mode b sin k x cos mt to u and 
the amplitude b is determined from the energy equation so that the. secular 
term associated with the mode (k,m) in u vanishes. Similar modifications should 
be made to avoid the appearance of a secular term of higher mode for the solu- 
tion of a forced oscillation of a nondispersive system. 

When the perturbation theory for the free oscillations yields only a 
trivial solution, we conclude that the small amplitude periodic solution which 
splits off from the state of rest does not exist. The longitudinal vibrations 
of a uniform bar with fixed or free ends are examples mentioned in ref. 1. This 
result is verified in ref. 1 by the method of characteristics . It yields a ^ 
stronger result that, for certain initial boundary value problems, every non- 
trivial solution becomes singular in a finite time. Shock waves are formed 
afterwards. This result is the same as that of Lax but the proof is somewhat 
different (ref. 6 1. This result was also applied (ref. 1) to show that there 
are no nonsingular (shockless) finite amplitude sound waves in a closed tube. 

For long tubes with open. ends, it can be extended to a problem periodic in x if 
the simple boundary condition of constant pressure is imposed. Consequently we 
can conclude ag^in the nonexistence of nonsingular periodic solution. However, 
the realistic boundary condition for ah open end of a long pipe should be im- 
posed based on the analysis of Levine and Schwinger (ref. 7). The reflection 
coefficient is not equal to -1. It depends on the wave number and furthermore its 
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absolute value is less than unity due to the loss of energy propagated to in- 
finity. Therefore an initial wave in an open tube will decay and a periodic 
free oscillation cannot be sustained. Detailed studies for this problem will 
be reported elsewhere. 
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